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Both spin and orbital degrees of freedom contribute to the magnetic moment of isolated atoms. 
However, when inserted in crystals, atomic orbital moments are quenched because of the lack of 
rotational symmetry that protects them when isolated. Thus, the dominant contribution to the 
magnetization of magnetic materials comes from electronic spin. Here we show that nanoislands 
of quantum spin Hall insulators can host robust orbital edge magnetism whenever their highest 
occupied Kramers doublet is singly occupied, upgrading the spin edge current into a charge current. 
The resulting orbital magnetization scales linearly with size, outweighing the spin contribution for 
islands of a few nm in size. This linear scaling is specific of the Dirac edge states and very different 
from Schrodinger electrons in quantum rings. Modelling Bi(lll) flakes, whose edge states have been 
recently observed, we show that orbital magnetization is robust with respect to disorder, thermal 
agitation, shape of the island and crystallographic direction of the edges, reflecting its topological 
protection. 


A central notion in magnetism is the fact that orbital 
moments associated to circulating currents are fragile. 
They naturally occur in open-shell isolated atoms [ij, but 
these atomic orbital moments quench as soon as the atom 
is placed in a crystal. Circulating currents in artificially 
patterned mesoscopic quantum rings 0, studied in the 
last three decades jj, 1^, require very special conditions 
to survive, such as very low temperatures so that the 
electrons keep their phase coherence around the entire 
ring, and small disorder, so that electrons do not local¬ 
ize. In contrast, robust spin currents occur naturally at 
the edge of quantum Spin Hall insulators (QSHI)[§-l3 
and are topologically protected. These spin currents are 
associated to Kramers doublets, where each state has a 
net charge current flowing with opposite chirality. In a 
finite sample, these counter-propagating currents can be 
associated to magnetic moments with opposite sign for 
each state in the Kramers doublet. Since these states 
are equally occupied, the resulting net orbital moment 
vanishes. Having an insulating bulk and robust spin cur¬ 
rents at the edges, QSHI are natural quantum ringsQ for 
spin currents. The central idea of this paper is that, in 
the case of QSHI nanoisladsQ (or flakes) with a discrete 
edge state spectrum, it is possible to turn these robust 
spin currents into robust charge currents that result in 
very large orbital moments. To do so, two conditions are 
sufficient: a magnetic field has to split the Kramers dou¬ 
blets and, using electrical gating or chemical doping, only 
one electron has to occupy the highest occupied Kramers 
pair, providing thereby a net edge current, and a large 
orbital magnetization. 


Several systems have been predicted to be QSHl[I3,[ni 
and strong experimental evidence exists that CdTe/HgTe 
quantum wells [l^, InSb/GaAs quantum wells [isj, and 


quai 

with Bismuth (111) monolayers |14l4l7j| host spin filtered 
edge states essential for our proposal. To substantiate 


our claim, we choose Bismuth for two reasons: the topo¬ 
logical edge states of nanoislands of Bi(lll) have been 
recently observed by means of Scanning Tunneling spec¬ 
troscopy (STM) [171, M, and a very well tested tight- 
binding Hamiltonian [l9| is available that makes it possi¬ 
ble to compute the electronic structure of systems with 
thousands of atoms. 


Bi(lll) bilayer (BL) is a buckled 2D honeycomb crys¬ 
tal (see figure [T|). We model the Bi(lll) nanostruc¬ 
tures with the same tight-binding model [l9| employed by 
Murakamito predict that Bi(lll) would be QSHI. 
The same approach has also been used by Drozdov et 
al. 0 and by Sabater et a/.[I^ to model their experi¬ 
mental results. The Liu-Allen tight-binding model de¬ 
scribes Bi with four orbitals {s,Px,Py,Pz) per atom, with 
interatomic hoppings up to third neighbors, parametrized 
with the Slater-Koster approach[2lj and atomic spin- 
orbit coupling XL ■ S. Within this model Bi(lll) bilay¬ 
ers naturally come as QSHl[i3, with spin-filtered edge 
states and a gapped bulk. The effect of the magnetic field 
is incorporated by using Peierls substitution 2^ , with an 
extra phase accumulated by electron going from site i to 
j, ipij = 27r-^ Adi, where A is the vector potential, 
^0 = — is magnetic flux quantum. 


The electronic structure of a Bi(lll) nanoisland, 
with hexagonal shape with six zigzag edges with length 
Ledge = 3.6nm each, is shown in Fig. [1^. In panel [T]d 
we show the energy levels corresponding to the 2D ma¬ 
terial, to mark the gap of 0.25eV, side by side with the 
discrete energy spectrum of the island. We denote the 
energy difference between adjacent Kramers doublets by 
5. Upon application of a magnetic field perpendicular 
to the island, the in-gap Kramers doublets split follow¬ 
ing straight lines, indicated by red (blue) colors for states 
with increasing (decreasing) energies. Red and blue color 
lines distinguish the in-gap edge states rotating clock- 




2 



FIG. 1: Orbital magnetization of Bi(lll) nanoisland edge states a Scheme of edge currents of a given 
Kramers doublet. At zero field, the occupation of both states with opposite orbital magnetization is the same (red 
and blue arrows). Application of a magnetic field, plus single occupancy of a Kramers doublet, results in net orbital 
magnetization (thick blue arrow), b Calculated energy spectra of a the 2D Bi(lll) bilayer and a flake with edge 
length L ~ 3.6nm (N=384 atoms) and evolution of flake spectrum as function of a magnetic field. The bilayer gap 
Egap, energy level spacing (5, and the splitting of the Kramers doublets in a magnetic field As are indicated. The 
corresponding magnetic moments, = | 1 are shown on the right. The largest magnetic moment M^ax is 

highlighted in green, c Calculated local current density (blue arrows) flowing along edges for the state generating 

Mmax • 


wise and counterclockwise, respectively. The splitting, 
denoted by As, exceeds by far the spin Zeeman split¬ 
ting, which actually is not included in the Hamiltonian. 
This indicates that these states carry an orbital moment. 
Its origin becomes apparent upon inspection of the plot 
of the current density associated to the one of two states 
of a Kramers doublet, displayed in figure (Tt, that shows 
the circulating edge current. It must be noted that some 
states not in the gap also have large orbital moments. 
We have verified that their wave functions are not fully 
localized at the edge, so that the emergence of the large 
orbital moment could have a different origin. 

To be more quantitative, we use the definition of mag¬ 
netic moment associated to a given quantum state [l| 


M„ = 


dEn 
dB ■ 


( 1 ) 


The absolute value of the corresponding orbital magnetic 
moment associated to the in-gap Kramers doublets is 
shown in mo. In Fig. [2] we show the magnetic moment 
for the in-gap state with the largest Mn for a given is¬ 
land, denoted as M„iax, as a function of the island size L. 


The magnitude of Mmax scales linearly with L, reaching 
values as high as 42 Bohr magnetons (/is) for a hexag¬ 
onal island with L = 18nm, much higher than the spin 
contribution (l/rs). For Schrodinger particles in a ring, 
the magnetic moment is given by Mschro = where 

Lz is the azimuthal quantum number, independent of the 
size. In contrast, from the exact solution of a massless 
Dirac particle moving in a ring 2^ 241, equation ([T|) gives 
that the magnetic moment of Dirac particles also scales 
linearly with the circle radius. Therefore, our results can 
be interpreted as if the edge states were described by 
Dirac particles confined in a ring. We can also relate 
the magnetic moment Mmax to a persistent edge current 
through the classical definition of magnetic moment in a 
loop Mniax = I A, where A is the area of the nanoisland. 
The current so calculated are shown in Fig. [5^, and can 
reach values of 0.5/iA for L = 18nm. 


Whereas quantum transport experiments could probe 
the magnetic moment of individual states [illii , mag- 
netometry experiments are sensitive to the total magne¬ 
tization, which involves contributions from all occupied 
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FIG. 2: Size scaling analysis, a Maximal orbital magnetic moments Mmax (see Fig. [5 }d) in units of Bohr magnetons 
(black squares) and corresponding current amplitudes Imax (red circles) as a function of edge length L of the flake. 
A linear dependence of orbital magnetic moments as a function of edge length L is clearly seen, with ~ 2.5^. 

b A comparison between characteristic energy scales in the systems as a function of edge length L. Energy level 
splitting of the edge states in the absence of the magnetic field 6 (blue circles) and a Kramers degeneracy splitting in 
a magnetic field Ag for a magnetic field B = IT (black squares) and B = 2T (red triangles). i5 = WimeV/L{nm) 
(blue curve) obtained from a fitting procedure, which is characteristic size-dependent quantization (oc L~^) for Dirac 

Fermions. 


states: 

( 2 ) 

n 

where the sum runs over all the state of the island and 
/(e„) is occupation of the individual states. In Fig. [3[i 
we show the T = 0 magnetization as a function of Fermi 
energy Ep of the island considered in Fig. [U for a mag¬ 
netic field oi B = IT. It is apparent that, whenever a 
split Kramers doublet is singly occupied (as seen in fig¬ 
ure WP: left panel), the net magnetization is very large 
and parallel to the applied magnetic field, corresponding 
to an orbital paramagnetic response of the island. When 
an extra electron is added or removed from this situa¬ 
tion (figure [3 )d, right panel), the total magnetization is 
small and antiparallel to the applied field, so that the 
island behaves diamagnetically. The total magnetization 
is roughly given by the magnetic moment of the highest 
singly occupied Kramers doublet. 

We now address the robustness of the orbital magne¬ 
tization with respect to thermal disorder, assuming that 
thermal equilibrium has been reached so that /(cn) = 
exp [ / 3 (£,!-/j) ] +i ’ where (3 = l/fcsT and ^ is the chemical 
potential, which we fix midway between the two states of 
the highest occupied split Kramers doublet. The evalua¬ 


tion of (Mtot) using eq. ([2]) requires the numerical calcu¬ 
lation of the entire spectrum, possible only for sufficiently 
small islands, such as those shown in Fig. 01 with L — 3.6 
and L = 4.5nm and (Mtot(T = 0)) = 5.8 and 7.5 /tb re¬ 
spectively. Upon heating, the magnetization remains sta¬ 
ble up to I Kelvin and then decays. The temperatures at 
which the magnetization decays by 50 percent are 3 and 
4 Kelvin respectively. Since the dominant contribution 
to the magnetization comes from the highest occupied 
Kramers doublet, with energies Ei,E 2 = Ei + As and 
magnetic moment M, we expect that the magnetization 
will be approximately given by: 

p ‘2kQ T _ p 2fc^ T 

(Mtot) « M {f(E,) - f{E2)) = M - 

2-|-g2fcBT _|_ g 2k^T 

The good qualitative agreement between the exact and 
the approximate curves, shown in Fig. 01 supports the 
use of the approximate equation to estimate {MtotiT)) 
for larger islands, for which numerical diagonalization is 
out of reach. Expectedly, larger islands have larger mag¬ 
netic moments, larger energy splittings Ab and therefore, 
the magnetization is more robust with respect to thermal 
occupation of states with opposite magnetization. For a 
Bismuth(lll) hexagonal island with L = I8nm our cal- 
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FIG. 3: Total magnetization, a Total magnetization Mtot as a function of Fermi energy Ep at temperature T = 0 
for the island with edge length L ~ 3.6nm considered in Fig. [TJ The off-gap states, as defined by the bilayer spectrum, 
are shown with a grey background. As the Fermi energy moves in-gap, total magnetization oscillates . b A schematic 
picture explaining magnetization oscillations shown in a. Red and blue bars correspond to energy levels from a given 
Kramers doublet split by a magnetic field. Filled states (yellow area) are occupied by electrons (indicated by arrows) 
in edge states. When number of electrons is odd, only one of the states from the highest Kramers doublet is occupied 
(left panel), inducing orbital magnetic moment +M that contributes to total magnetization. In this state, addition 
or removal of a single electron results in a quenching of the edge magnetization. 


culations predict an orbital magnetization at T = 4K as 
large of 42^^, that would only be depleted by 10 percent 
at T = lOK. 

Two energy scales determine whether is possible to 
selectively occupy a single state in a Kramers doublet: 
the energy splitting between different Kramers doublets 
(5 and, within a given Kramers doublet, the magnetic 
splitting As. As we show in Fig. [^b, 5 oc which 
again reflects both the edge character of these states as 
well their Dirac nature, whereas oc L, as expected 
from the linear scaling of with L, shown in Fig. [2^. 
So, increasing the size of the islands makes the magnetic 
moment of individual edge states larger, but eventually 
makes it impossible to prevent scrambling of Kramers 
doublets. Therefore, there is a magnetic field dependent 
optimal size 6{L) ~ Ab{L) for which orbital magnetiza¬ 
tion is maximal. We note that maximum orbital moment 
as a function of applied magnetic field is stable as long 
as there is no crossing with other states, at values of B 
so high than the magnetic splitting As is larger than the 
zero field splitting S (see Supplementary Materials). In 


that case, two states of different Kramers doublets could 
anticross resulting in a drop of their orbital moment. 

The robustness of edge spin currents in QSHI is due, 
ultimately, to time reversal symmetry Q. Time reversal 
symmetric perturbations can not produce elastic edge 
backscattering. Application of a magnetic field breaks 
time reversal symmetry, which combined with a time re¬ 
versal symmetric disorder potential could, in principle, 
produce backscattering, resulting in the destruction of 
the orbital magnetization by mixing states with opposite 
orbital magnetization. Therefore, we test the robustness 
of our predictions by studying the effect of disorder. We 
first consider Anderson disorder, introducing a uniformly 
distributed random potential on every orbital of the sys¬ 
tem. This introduces both atom to atom variations, the 
conventional Anderson disorder, but also random crystal 
field splittings at every atom. 

The average T = 0 total magnetization for B = IT 
for the island considered in Figs. [T] and [31 obtained after 
after averaging over 100 realizations of disorder config¬ 
urations, with on-site energies randomly distributed on 
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FIG. 4: Temperature dependence of total magnetization. Temperature dependence of total magnetization 
Mtot for systems with different sizes, L ~ 3.6nm (green lines), L ~ 4.5nm (red lines) and L ~ 18nm (blue line), for a 
magnetic field B — IT and chemical potential in the middle of two states from Kramers doublet (see left panel in Fig. 
[3b). Solid lines correspond to results obtained using eq. m and dash lines to two-level approximate model given by eq. 

For zero temperature the state with orbital magnetic moment +Mn contributes to total magnetization M. With 
increasing temperature, the state with opposite orbital magnetic moment —is populated, reducing magnetization. 
The blue dash line predicts magnetization of the system with edge length L ~ 18 nm by using eq. ©■ No change 
of magnetization up to T = 5K is predicted which is related to large energy level separation Ab{B = IT) r; 5meV 

shown in Fig. EKb). 


the interval ±bF/2, denoted by {{Mtot)){T = 0) is shown 
in Fig. ini We also plot the statistical standard devia¬ 
tion, but due to its small size it is smaller than the data 
points. The stability of the orbital magnetization is re¬ 
markable even for disorder strength of W = 0.5 eV per 
atomic orbital. The very week effect of disorder on the 
edge magnetization can also be seen (inset of Fig. |S]) 
in the evolution of the in-gap edge states spectrum as a 
function of the applied held B. 

In addition, we have verihed that positional disorder at 
the edges (see supporting information Fig. Sic) does not 
reduce the orbital magnetization either. We have also 
verihed that the shape of the island and the type of edge 
play no role: similar results are obtained for triangular 
zigzag islands and for hexagonal armchair islands (see 
supporting information). 

We have also considered the inhuence of the substrate. 


relevant for the case of Bi(lll) hakes on top of Bi(lll) 
bulk[I3. For that matter we have calculated the orbital 
magnetization of the edge states of zigzag hexagonal is¬ 
lands on top of a much larger Bi(lll) hake. Our results 
show that the orbital magnetization of many edge states 
is still preserved in the supported islands(see supporting 
information). 

The phenomenon of robust orbital magnetization in 
nanoislands of QSHI can be also analyzed from a dif¬ 
ferent perspective: these systems behave like mesoscopic 
quantum rings fabricated to observe persistent currents. 
These QSHI nanoislands have three major advantages, 
compared with the conventional quantum rings. First, 
there is no need to pattern any structure, since the bulk 
of QSHI islands is not conducting. Second, the topolog¬ 
ical protection of the edge current how, which our cal¬ 
culations show is preserved in part in the presence of 















6 





> “I-'-1-'-1-'-1-'-1-'- 

0.0 0.1 0.2 0.3 0.4 0.5 


W[eV] 


FIG. 5: Robustness of magnetization with respect disorder. Stability of magnetization Mtot against disorder 
W at temperature T — 0 for a nanoisland with edge length L ~ 3.6 nm. A random onsite energy per atomic orbital a 
is chosen from an energy interval Ea = (—VF/2, Wl2). Magnetization M decreases by less than 10% for the strongest 
disorder, W = 0.5. The plot includes also statistical standard deviation bars, however due to their small values they 
overlap with data points. The inset shows evolution of energy spectrum in a magnetic field B for disordered system 
with W = 0.2 eV. Dispersion of edge states in a magnetic field B is not significantly different in comparison to clean 

system, shown in Fig. [T] 


magnetic helds and disorder, results in a robust persis¬ 
tent current. Third, the Dirac nature of the quasiparti¬ 
cles permits to upscale the resulting magnetic moment 
linearly with size, unlike Schrodinger quasiparticles, for 
which the magnetic moment does not depend on size. Ac¬ 
tually, the existence of edge states with persistent charge 
current was also predicted for QSHI quantum dot made 
of HgTe quantum well with inverted band structures 271. 
We have also verified (see Fig. S3 in Supp. Mat.) that 
robust orbital nanomagnetism is present in the Kane- 
Mele model for a nanoisland. Our results, together with 
those of reference [13, suggest that the physics discussed 
here is universal, model independent, and thereby can be 
expected from any QSHI nanoisland. 


Experimental work will determine if the orbital mag¬ 
netization will present remanence, as in the case of nano¬ 
magnets, or at zero held the orbital moment will present 
random telegraph noise as in super paramagnetic par¬ 


ticles. This will be related to the fascinating question 
of spin relaxation between two states in a Kramers dou¬ 
blet, that entails a rather large change in orbital angu¬ 
lar momentum[28|, and coupling to other spin degrees 
of freedom, such as the Bi nuclear spins, will certainly 
play a role. The selective occupation of a spin havor at 
B = 0, and the resulting orbital magnetization, could 
also be driven by Coulomb interactions [2^ that we have 
neglected in this work. Finally, it has not escaped our 
attention that these islands could operate as well as very 
good interface for single spin readout: injection of a sin¬ 
gle electron in an otherwise closed shell conhguration will 
result in a orbital magnetization conditioned to the spin 
orientation of the added electron. The magnitude of the 
orbital magnetization would be well within the reach of 
state of the art local probe for magnetization, such as 
magnetic resonance force microscopyjs^ and NV center 
nanomagnetometry . 
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The edge type and nanoisland shape dependence 

We have verified that the orbital magnetization occurs regardless of the edge type and nanoisland 
shape. In Fig. [^we show an evolution of energy spectra in a magnetic field of a triangular nanois¬ 
land with zigzag edges, b hexagonal nanoisland with armchair edges, c hexagonal nanoisland with 
structural edge disorder. All these spectra are calculated using the same four-orbital tight-binding 
model discussed in the text.- Whereas the distribution of level spacings d of the islands with arm¬ 
chair edge and the one with a disordered edge is less regular than in the case of islands with ideal 
zigzag edges, all of them show a strong orbital magnetization, as reflected by the large splitting 
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induced by a magnetic field. Results presented in Fig. [^are in agreement with an analysis from 
Ref.^ regarding an existence of topologically protected edge states for structures with arbitrary 
edge orientation. For the island with dangling atoms a linear splitting in a magnetic field for many 
in-gap states is still observed. These calculations confirm the presence of orbital magnetization in 
arbitrary shape and type of edges QSHI nanoislands. 

a b c 




Figure S1: Shape effect. An evolution of energy spectra in a magnetic field of a triangular nanois¬ 
land with zigzag edges, b hexagonal nanoisland with armchair edges, and c hexagonal nanoisland 
with structural edge disorder. For all structures some set of states reveal linear splitting in a mag¬ 
netic field. 

Influence of the substrate 

We now address the influence of the interaction with the substrate on the orbital magnetization. 
Motivated by the experiments by Drozdov and coworkers,- where the Bi(l 11) flakes were located 
on the surface of a thin film of Bi(l 11), we compute the spectrum for an hexagonal nanoisland such 
as the one in the main text, with L ~ 3.6nm {N = 384 atoms), on top of a much larger Bi(l 11) island 
with dimensions 50 x 55nm {N = 31240 atoms). Periodic boundary conditions are assumed along 
one direction (see Fig. [S2h l so that the substrate edges are disconnected. Full diagonalization 
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of Hamiltonian matrix of such a big system (in this case size of Hilbert space ~ 3 * 10^) is not 
possible, so we find only a set of eigenstates in a vicinity of bilayer energy gap Egap using iterative 
eigensolver method.- The energy splittings for B =\T, both with (inset) and without (main panel) 
coupling to the substrate, are shown in Fig. [S2b . We focus on the energy range given by the in¬ 
gap region, taking the gap Egap, from the two dimensional Bi (111) bilayer, which is marked by 
vertical red dash lines. We find additional in-gap states of the supported island, compared to the 
freestanding one, that correspond to substrate edge states, as inferred from inspection of their wave 
functions. In the inset we assign a color code to the states: red stands for substrate states, blue for 
island states. Whereas the evolution of the magnetic moments as a function of energy is no longer 
a smooth function, it is still apparent that many in-gap states retain their orbital magnetization. The 
substrate in-gap states have a very small splitting (red points). Thus, our calculations suggest that 
orbital edge currents could still be present in Bi(lll) flakes deposited on Bi(lll). We have also 
verified that similar results hold for larger islands, with L ~ 9 nm edge length {N = 2400 atoms) 
deposited on the same substrate. 

A magnetic field dependence 

A magnitude of maximum orbital moment as a function of applied magnetic field is investigated in 
Fig. [S3l The inset shows a corresponding energy levels evolution in a magnetic field (a green line 
indicates a state generating Mmax)- The magnetic moment is stable, only slightly decreasing, up to 
B = 70T. When a crossing of a given state with other state occurs, e.g. for B = 40T, see the inset, 
there is no change of orbital moment due to no change of a slope of energy level in B. On the other 
hand, when two states anticross, e.g. for B = SOT, it results in a drop of the orbital moment. When 
anticrossing region is left, magnetic moment return to its high value. 

Orbital nanomagnets with other Hamiltonians 

In order to check universality of orbital nanomagnetism in QSHI nanoislands, we model them with 
the Kane-Mele Hamiltonian.^ This model, widely used to describe the QSH phase, is a good ap- 
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Figure S2: Interaction with substrate, a Nanoisland with L ~ 3.6nm {N = 384 atoms) deposited 
on 50 X 55nm {N = 31240 atoms) substrate from top(upper panel) and side(lower panel) view. In 
order to reduee edge effeets from the substrate, periodie boundary eonditions in one direetion were 
applied. Energy splitting of Kramer’s degenerate pairs for a magnetie field 5 = IT without 
the substrate is shown in a main panel and with the substrate in the inset. The Bi (111) bilayer 
energy gap Egap is indieated by red dash lines. The states loealized within nanoisland are found by 
ealeulating eleetronie probability densities for eaeh energy eigenstate. States fully loealized in the 
substrate are indieated by red eolor while states fully loealized in nanoisland by blue eolor. 
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Figure S3: Magnetic field dependence, a Maximum orbital moment Mmax as a function of applied 
magnetic field for nanoisland with L ~ 3.6nm (A^ = 384 atoms). The inset shows an evolution of 
energy levels as a function of a magnetic field. A green line indicates an energy level generating 

^max- 


5 
















proximation for materials such as graphene and Silicene.- It features a single orbital per atom and 
the spin-orbit interaction is given by a spin-dependent second neighbor hopping. In the presence 
of an external magnetic field the Hamiltonian is written as 


HKM = t Y, e“^‘ja]^aja + i^so Y ( 1 ) 

with summation in the first term over nearest neighbors and in the second term over next nearest 
neighbors, and v,j = ±1 for clockwise and counterclockwise direction of path connecting site i 
and j, is a Pauli matrix and A 50 strength of effective spin-orbit coupling. We have simulated a 
hexagonal nanoisland with armchair edges consisting ofN= 1986 atomic sites and L = 4.4nm. In 
Fig. [S4]the evolution of the energy spectrum as a function of magnetic field is shown. Again, our 
calculations yield a linear splitting of Kramers doublets with a magnetic field around energy E = 0. 
Using the same argumentation of the main text, this should lead to the generation of persistent edge 
currents whenever a single electron occupies the highest occupied Kramers doublet. In the inset, 
we show the size-scaling of maximal magnetic moment M^ax from these states, and we find a 
linear scaling, consistent with the behavior expected for Dirac electrons in a ring. 

Finally, we note that quantum dots described with the Hamiltonian proposed by Bemevig- 
Hughes-Zhang (BHZ) to predict the QSHI phase CdTe/HgTe quantum wells,- as well as in III-V 
type II structures,- also lead to the existence of in-gap edge states^ in the dots endowed with 
orbital magnetization. Therefore, the concept of orbital nanomagnet built with zero dimensional 
nanostructures of quantum Hall insulators is a model-independent prediction. 
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Figure S4: Nanoisland within Kane-Mele model. An evolution of energy spectra in a magnetic 
field of hexagonal nanoisland with armchair edges with N = 1986 atoms and L = 4.4nm calculated 
within Kane-Mele model for t — — 3.0eV and Xso = lOmeV. In a vicinity of the middle of the en¬ 
ergy spectrum, £ = 0 a set of linearly dispersed states in a magnetic field is clearly seen. The inset 
shows a linear dependence of maximal orbital magnetic moment My^ax as a function of nanoisland 
edge length L. 
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